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ABSTRACT 


In  the  first  part  of  this  work,  the  isobaric-isothermal 
loop  of  the  methane-carbon  dioxide-hydrogen  sulfide  system 
at  the  pressure  of  1330  PSIA  and  temperature  of  160°F.  was 
investigated,  in  an  effort  to  complete  the  work  done  on  the 
phase  behavior  of  the  above  system  at  the  mentioned 
temperature. 

The  experimental  procedure  consisted  of  preparing  a 
known  mass  of  gas  mixture,  of  definite  composition,  with 
subsequent  compression  of  the  system  to  the  two  phase 
region.  Experimental  phase  composition  determinations  were 
carried  out  defining  thus,  the  three  component  two  phase 
region . 

The  experimental  data  were  correlated  by  means  of 
calculation  of  equilibrium  ratios  for  the  pure  components, 
and  the  deviation  from  ideality  for  the  liquid  and  gas 
phases  was  illustrated  by  means  of  activity  coefficient 
calculations . 

Comparison  of  the  experimental  results  obtained  to  those 
calculated  through  the  use  of  experimental  data  on  the  same 
system  at  some  other  pressures  and  temperature  of  160°F., 
indicates  substantial  differences. 

In  the  second  part  of  this  work  the  volumetric  behavior 
of  several  methane-carbon  dioxide-hydrogen  sulfide  mixtures 
was  investigated  at  temperatures  of  160°,  100°  and  40°F. 
and  pressures  up  to  3000  PSIA.  Three  carbon  dioxide-hydrogen 
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sulfide  binary  mixtures  were  also  investigated  at  the 
temperature  of  160°F.  and  pressures  up  to  approximately 
2000  PSIA. 

The  technique  employed  consisted  of  preparing,  in  the 
gas  phase,  a  mixture  of  known  composition  with  subsequent 
compression  and  recording  of  the  variables,  i.e.,  pressure, 
temperature  and  volume.  Correlation  of  the  experimental 
data  was  achieved  by  calculating  compressibility  factors 
for  the  mixtures  at  each  temperature  as  well  as  residual 
volumes. 

The  experimental  results  indicate  that  for  the  higher 
temperatures  the  compressibility  of  the  mixtures  can  be 
relatively  accurately  described  by  means  of  the  pure 
components  using  methods  reported  in  the  literature.  Some 
of  the  prediction  methods  fail  to  describe  the  compressibility 
even  at  the  higher  temperatures.  For  the  lowest  temperature 
of  40®F.  where  even  at  low  pressures  carbon  dioxide  and 
hydrogen  sulfide  are  in  their  two  phase  region  there  is  no 
adequate  method  of  describing  the  system's  behavior. 

It  should  also  be  mentioned  that  greatest  discrepancies 
occur  for  mixtures  poor  in  methane  content. 
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the  reciprocal  of  the  slope  of  the  isometric 
of  a  gas  as  defined  in  the  text. 

the  internal  pressure  of  a  gas  as  defined  in 
the  text . 

mole  fraction  of  a  component  in  the  liquid  phase, 
mole  fraction  of  a  component  in  the  vapor  phase, 
compressibility  factor  in  the  gas  phase, 
residual  volume. 

Subscripts 
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referring 

to 

liquid  phase. 
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referring 

to 

vapor  phase. 
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any  component. 
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GENERAL  INTRODUCTION 


The  past  twenty-five  years  in  petroleum  production 
and  processing  have  seen  a  broadening  in  the  range  of 
operating  temperatures  and  pressures.  To  produce  and 
process  hydrocarbon  mixtures  efficiently  at  these  varied 
conditions  has  required  the  extensive  efforts  of  many 
researchers  to  accumulate  and  correlate  basic  vapor-liquid 
equilibrium  data. 

Similarly,  all  petroleum  engineering  problems  use 
the  PVT  behavior  of  fluids  in  one  form  or  another.  Gas 
reserves,  gas  pressure  gradients  in  well  bores,  and  the 
calculation  of  gas  flow  rates  at  the  surface  are  three 
example  problems  requiring  a  knowledge  of  gas  behavior. 

An  accurate  knowledge  of  the  fluid  density  is  also 
required  in  the  calculation  of  orifice  meter  correction 
factors  in  the  metering  of  high  pressure  gases,  the  times 
of  contact  of  vapors  in  a  cracking  furnace,  and  the  vapor 
velocities  in  bubble  towers. 

Finally,  the  tendency  for  condensation  of  some  non¬ 
hydrocarbon  compounds  such  as  carbon  dioxide  and  hydrogen 
sulfide  -  present  in  a  variety  of  natural  hydrocarbon 
mixtures  -  which  is  favored  by  the  conditions  during  the 
refining  operations,  necessitates  an  adequate  knowledge  of 
phase  and  volumetric  behavior  for  these  mixtures. 


' ■'S' 

;  .  ■! 

.a  ,...;;:  ^'.rsi  uv-lc-AJ-  .vJ-  oTi's.v  .v/.i'-vdn:r;/;;  J-a^c:  ‘v.iT  '-V 

.-v:  .j.-,-.  .a/  ,,..1......  ]:  ..  ■■  i:  -uac'a  T.uixaaar->v<:::q  ^na 

,  .■  -.7; '■■.■•a  cq/  ,  .■,a',..;fiat, . ■  Lna.anrz/j'  .cuqo 

l7..L-.';'.-v  -.a.  -  '  ,.  ..A7l:£:,:'DO(:vqqu  ar-a: :..;o:r;q 

,.;.  1.7'  .a’r'"':a'7  avJ.c  i^caXv;.  ;:'q.  ;:y-iva::j^:  '  ::  r/i..  ^£7 ’.crq".  a^,.'D 

a;-,.  .L.'l 'M-m-:  ■ L  xa..-^'::  -/a a.*::  ;  vcxx'a  *  ■£!!:  .qy...Xi.:..L>  .  ‘J.j  axsTtax-^OBDi: 


ax.ar  a.:. 

: .  .da..’ '.a:'  xxa  a..3x  a.':. 

X  -n  ■ 

^  .Xc  ad  a  a  .. . ,  :  ,  a-x  a  a  1  .i:  raX  3  ' 

aa...  .a 

.:..L-..^ X  xx.'-.  :r.L  ■  la’i!':.:  ■ 

a  no 

nl 

i7fd;-x.:[?.  1.0  'icx'/x'xi ..  d  .'  a.. 

ra'iJ- 

lPP  I,£..-..  ,::::>v;C''  J.I: -w  nx  a .  x.-.q?  ..■'v;ic  aaKX  b.',;  ,  ar  v  7;;:>a  .av 

.:■::  7;.a:J'  .xa  a..-  V;v,.jr,  ,.u<'.j  JiS  a  iaa.a:  o/>p  “•'•■'  X‘).,.:-i'i..-.-.i:r;:>.L7.io 


.  xcx  /..ii..id'  r,xp  :.(■.•  xpbc-: X'wf' x:;: .'  .a  ■.piixa.lupa x  a(;icIc[o'i«r  '.;.jP|in.xxo' 
r  ;;.  J.  ,1  ax  'v  aa.a’.i.aa'jr;  £.:a’^£_  ~-.^..is  x.*  .  a q,: .x; X iiX  ■:■.£  -'  V-jajaa.  u,. ; 

.a:  .'aj-.j;-'’.. '.-vi.  a:,.:;.-f!T  DQi;:d:aa>  Xo  xadda'aCB...  .a...  an:,  ad;i.ii;pa'i 

::.a,,ij:a  .:x[a  ,  - ^-'a-  a-j.^r';:,aaxa  iI’Diff  ao  pa;.:aadBta  Bria  xJx  aa-odoXa 

'.a  <  ij  V  .  ijy  .  :..  ^  .  ..  xaaaa.  .jx ..pl'a a.aa  a  d:.!  axocpaX  ao  a  xico  Xo 

»jaawc-X  .'Ii'Jja:  ad:  aaatdXaca^ev 
".rax.  i-x  X.,  '‘X  iia'.d/ •  :;aa.-j' xa  ^1.'  'xc^  v^aBaai^ea;  a:>£x.  ^  pAXiaaid; 

..a;n  .' :  V  i  .;' j  ..■_  ..  i'  a...  t  .':i'  i::o  :zi..:.<  i:c:  xxxaa-a  cj'; j.ii rqqfacQ  .a* '.la:.....  c 

ii-  dx:.ao:u..va;  .;;  v.. x...  .  xx  :  ■  ■  q:  .'J;a  :v  .  xa'-.  ...  -  'a:..-.'::  .d/a 

'-..a.  ^■x’i.ia;ai  ..  j..  .L/J:'.  Xi'a.,  ..  j.:,  '/rX  a  a-  a.£  ^xx.^iilv/-  •-.  .i,0':a/aa:j?m 

j  ^.r.j  j.  '  : .  .  a.  :'f  ,.■  .  ki.:>  c.;  Jaa' j.:.a-.:  aaa.;  \v.ixc:'.:.a.  'laoo  x.'.aaxi,^  ^  sa 

...  ...ixJ- i:.a  ■  .i.:.xxJ  axa:  axahv  ■xOa,  :;:.±:i:a7.;f  aaf^av  CaX..;  aa  '..dq 


-  ii 


Although  phase  equilibriuin  and  compressibility  factor 
measurements  of  numerous  binary  systems,  containing  a 
paraffin  hydrocarbon  and  carbon  dioxide  or  hydrogen  sulfide, 
have  been  published,  there  is  only  a  limited  number  of 
papers  dealing  with  the  ternary  system. 

The  work  of  J.  A.  Bailey  (1)  and  A.  P.  LorenZo  (2) 
at  the  University  of  Alberta  has  contributed  to  the  phase 
equilibrium  study  of  the  methane-carbon  dioxide-hydrogen 
sulfide  system,  at  three  temperatures  and  several  pressures 
up  to  1800  PSIA. 

The  aim  of  the  present  work  is  to  complete  the 

investigation  of  the  above  ternary  system  by  determining 

its  isobaric-isothermal  section  at  1330  PSIA  and  160°F.^ 

and  to  study  the  volumetric  behavior  of  the  system  in  the 
(vapor) 

single^phase  region  at  various  compositions,  temperatures 
of  40,  100  and  160°F.  and  pressures  up  to  3000  PSIA. 
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PART  I.  PHASE  BEHAVIOR  OF  THE  METHANE- 


CARBON  DIOXIDE-HYDROGEN  SULFIDE 
SYSTEM  AT  160°F.  AND  1330  PSIA. 


A. 


INTRODUCTION 


The  primary  purpose  of  the  first  part  of  the 
investigation  was  to  extend  the  study  of  the  methane- 
carbon  dioxide-hydrogen  sulfide  phase  behavior  relation¬ 
ship  begun  by  J.  A  Bailey  and  A.  P.  Lorenzo;. 

It  was  felt  that  for  the  completion  of  the  system 
at  the  temperature  of  160®F.,  an  additional  isobaric 
section  at  a  pressure  in  the  neighborhood  of  1300  PSIA 
should  be  studied.  The  reason  was  an  uncertainty  as  to 
the  shape  of  the  K  vs  P  curves  for  each  of  the  three 
components  at  the  mentioned  pressure  region.  It  should 
be  mentioned  here  that  A.  Lorenzo  covered  the  isothermal- 
isobaric  section  at  160®F.  and  1300  PSIA  by  employing 
prediction  methods,  utilizing  data  at  some  other  pressures 
and  the  temperature  of  160®F. 

No  modifications  of  the  existing  apparatus  were  made 
and  the  same  general  experimental  technique  and  sampling 
methods,  followed  by  Bailey  and  Lorenzo >,  were  used. 

The  data  were  analyzed  and  compared  to  the  figures 
calculated  by  Lorenzo^  for  the  predicted  isothermal-isobaric 
section  of  160®F.  temperature  and  pressure  of  1300  PSIA. 
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B.  SOME  THEORETICAL  ASPECTS  ON  THE  PHASE  BEHAVIOR 

OF  A  MULTICOMPONENT  SYSTEM 


1.  Historical  Consideration 

Early  investigators  discovered  that  at  relatively 
low  pressures,  mixtures  obeyed  relatively  simple  rules, 
e.g.,  Dalton's  and  Raoult's  laws  combined:  K  =  vapor 
pressure  divided  by  total  pressure,  and  therefore 
prediction  of  vaporization  ratios  for  hydrocarbons  was 
generally  believed  to  be  a  simple  matter,  i.e.,  = 

Yi/Xi- 

Later  investigators  discovered  that  at  high  pressures 
additional  considerations  were  required  for  predicting 
vaporization  ratios,  e.g.,  K's  defined  in  terms  of  pure 
component  fugacities  which  gave  good  results  up  to  several 
hundred  PSIA.  Both  types  of  K's  were,  by  definition, 
independent  of  composition.  That  is,  for  a  given  component^ 
K  the  equilibrium  ratio  y/x  was  defined  as  a  function  of 
temperature  and  pressure  only. 

Still  later  investigators  learned  how  these  simple 
rules  actually  break  down  progressively  with  increasing 
pressure^  and  become  totally  inadequate  at  pressures 
approaching  the  critical  pressure.  At  these  pressures  K 
values  apparently  depend  to  a  large  degree  on  composition 
as  well  as  temperature  and  pressure;  that  is,  the  K's  for 
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any  component  depend  also  upon  the  presence  and  the  amounts 
of  the  other  constituents  in  the  mixture. 

Also,  pronounced  effects  on  K-values  were  noted  in 
systems  containing  compounds  differing  from  the  paraffins 
in  molecular  structure,  such  as  aromatic,  as  well  as  non¬ 
hydrocarbons,  such  as  alcohols,  H2S,  CO2 /  N2/  etc.  In  short, 
K-values  are  not  constant  but  vary  for  each  component  depending 
on  temperature,  pressure  and  composition  of  the  mixture. 

To  obtain  quantitative  inf ormation  related  to  the 
effects  of  composition  on  K-values,  numerous  studies  v/ere 
made  covering  the  entire  two  phase  regions  of  various  binary 
hydrocarbon  systems,  of  several  ternary  mixtures  and  some 
more  complex  mixtures.  In  each  case,  a  pronounced  effect  on 
K-values  was  measured  as  the  mixture  approached  its  critical 
conditions.  Many  attempts  were  then  made  to  define  quantitatively 
this  radical  departure  from  ideal  liquid  and  vapor  behavior. 

One  approach  was  to  define  vapor-liquid  behavior  in  terms  of 
thermodynamic  relationships  and  PVT  behavior  of  pure  components. 
Such  attempts  are  the  arbitrarily  derived  equations  of  Van  Laar 
and  Margules. 

This  approach  which  has  promise  of  a  high  degree  of  accuracy 
is  impractical  for  routine  calculations.  A  second  approach  was 
to  utilize  the  large  volume  of  directly  measured  vapor-liquid 
equilibrium  ratios  which  has  accumulated.  This  approach 
makes  use  of  the  concept  of  convergence  pressure  to  describe 
the  influence  of  composition  on  equilibrium  ratios. 
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However  the  greater  part  of  this  information  has  been 
collected  for  binary  hydrocarbon  systems,  whereas  the 
primary  need  has  been,  and  is,  for  vaporization  ratios 
of  hydrocarbons  in  multicomponent  systems  and  mixtures 
containing  non-hydrocarbon  compounds.  Therefore,  although 
an  extensive  experimental  program  for  a  multicomponent 
system  is  time-consuming^  it  is  the  best  means  now  available 
for  the  determination  of  phase  equilibrium. 


2.  Three  Component  Phase  Equilibrium 

The  methods  of  representing  the  ternary  system 
graphically  are  sufficiently  different  to  warrant  more 
explanation  than  is  necessary  for  the  binary  system.  The 
phase  rule  statement  of  the  equilibria  which  may  exist 
in  a  ternary  system  is  summarized  in  the  following  formula: 

F=3+2-^  =5-0  where  F  the  degrees  of 

freedom  of  the  system  and  0  the  number  of  phases. 

For  a  single  phase  system,  when  F  =  4,  i.e., 
temperature,  pressure  and  two  compositions,  various 
methods  have  been  proposed  to  represent  the  composition  of 
ternary-isobaric  systems  which  are  based  on  the  triangle, 
but  only  two  have  been  used  to  appreciable  extent. 
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Since  phase  equilibrium  measurements  are  mostly 
carried  out  under  isothermal  conditions,  the  temperature 
is  selected  as  the  invariant  quantity.  Then  a  three 
dimensional  representation  of  the  phase  equilibrium  can  be 
acoomplished  by  selecting  pressure  as  the  ordinate  and 
two  compositions  as  abcissas.  The  so-constructed 
three  dimensional  model  can  be  studied  as  a  number  of 
parallel  planes  -  isobaric  sections,  of  the  form  of  an 
equilateral  triangle. 

The  first  method  proposed  by  Gibbs  (3) ,  makes  use 
of  an  equilateral  triangle  of  unit  altitude.  The 
quantities  of  the  components  are  expressed  as  fractional 
parts  of  the  whole;  the  sum  of  the  fractional  parts  is 
therefore  equal  to  unity. 

The  corners  of  the  triangle  represent  the  pure  components. 
A  point  on  the  side  of  the  triangle  represents  the  composition 
of  a  binary  system,  and  a  point  within  the  triangle  represents 
a  composition  of  the  ternary  system. 

The  second  method,  due  to  Roozeboom  (3 ) ,  has  found 
greater  usage.  An  equilateral  triangle  is  used  but  the 
length  of  the  side  is  made  equal  to  unity.  The  composition 
of  a  ternary  system  is  determined  by  this  method,  by  the 
distances  to  the  sides  in  a  direction  parallel  to  the  sides, 
in  distinction  from  the  first  method  where  the  perpendiculars 
from  the  vertices  to  the  sides  are  used.  Some  useful  features 
of  the  triangle,  whether  equilateral  or  not,  should  be  noted: 
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1.  The  ratio  of  two  components  is  constant  along  a 
line  through  the  third.  The  composition  of  a 
point  then  may  be  expressed  either  in  terms  of 
%A  and  %B  (any  two  components)  or  in  terms  of  °/C 
and  the  ratio  a/b  as  the  independent  composition 
variables . 

2.  The  ratio  of  A  to  any  fixed  composition  such  as 

D  on  the  opposite  side  is  similarly  constant  along 
a  line  through  either  of  the  other  corners.  This 
principle  can  be  extended  whenever  necessary  or 
convenient  for  a  sub-triangle  in  the  actual  ABC 
triangle. 

3.  If  two  complexes  or  compositions  are  mixed,  the  total 
composition  lies  on  the  line  joining  them,  at  a 
point  dividing  that  line  into  segments  inversely 
proportional  to  the  ratio  of  the  quantities  of 

the  original  complexes  taken. 

Some  of  the  above  properties  are  used  in  a  graphical 
construction  of  the  equilibrium  loop  from  experimental  data 
relative  to  an  arbitrarily  chosen  composition  parameter  C, 
as  is  explained  in  Lorenzo  ‘s  and  Bailey's  works. 

It  is  worthy  of  note  that  each  side  of  the  as  above 
formed  isothermal  prism,  represents  the  pressure  composition 
loop  for  one  of  the  binaries?  that  bubble  and  dew  sections 
are  represented  by  two  different  surfaces  inside  the  prism? 
that  the  composition  of  any  two  coexisting  phases  are  represented 


t. 

fi:  ■ 


r.-:  :-iii>-:J-c.nr,:::)  £:;i:  :  :  x....iAi;cciOo  c,.v/:;  -.c  p 

/.  ;  .  ’;iXo^  qr  ;;  vrf:,.  ..„'.L:uq  rfP0CxifJ 

;  ...  Pi.iuPi  ii::;  •f'OE-xx’ :i<:p:p  .='<  'V.effi  :;i-.EXq 

^  -...!:  ■-;■  ^  (...■.:l.i'i:: -r-oqiViOr.  cv,::;  '  ■  -•■  ■ 

i;.  x  aV.u:u.1:xx':  -ax  t;\A  c\-rx.x  .-dj 

.1- I.  :xj:-.  ■/ 

:-:  riL^^x...,  ....:  .:..-xx  ^;  v;.r:^  c; :.:  ^  :x.-'  .x^-.  .;  --i'.'- 

■j.,.  j  r  .At;.  ..  -  ^^x;:.I.i:ii:a.^,  x..x  x.x..[:o'  .  dd  ac  a 

.•,-;.,r;q;o  "jLrLijo  Bid?  z.o  X-  lipiXv: 'i.r:;  rxnxi  B 

r  ,  ., 


jr.  ■  n  X;.-' V::  X'-  ivV'  .OXX.' uiX  X Ix-S  OX.  O  ^  X -£4 

L>c'£/  S.  ^y  -:  -rIJ-  n.i:  o- 1 p no xxO a  B  oo±  ■ oOj.novxi' o> 

'■  ■  „  oJ[  "'OXT 

rxJ'OD'  oflx  ^  .'Ixo .  .Oil  OX'';  c  ix. lax -qr/iOo  oo  a... .'vO..!  p;  :c.o  -  .x  -o 

■J3  00  /.o.' X..U.  r.  xo.a  xL:  .00  o.  ,  .X  ..o-,q^Ox..x  ._;xo 

xx’; .' ax;:.' Viil  ol'.uOiTrOwO  cx-.ul  .  n/.!  drddd  'Oiix;  ■xvx.b  x.oxv. a 

:.a  c.  caxr.oo]:  ^  rib  xc  <.  .:  o .>•.:,  aril  OO'  I.xr'^  .ix  o^'X.:-.oq  " 

,  .aixx'b  0X0 1.4000  J.,■o.;.;.p.i: oul'  ■ 
p  ',.  .--.pv  x,b.  ooaa  co-O  oe.xxxpc[C  a; -x  r  vooo  rrib  icx  oror  L 
.ab.  b  Jo  Oil!  ■i.;l'0..(o:,o  0100:1:  cro  I  ■  :xi7.toc,'.i:X.oipa  or^b.'  :o.  1.0,  i:.iD;!:i':cbc,ii(a>o 


■.t  1  i/:o o:i  loor-  r;.o  X’  .:..  vli x  Jo.  o 


■OX'. 


1:  bns 


;  vl- :  L 

..K...;.;:-.'.  aoc:;.:o  r 


.......i,:oo(.  t.I  .a..X;  jD^xoi:.:  IX'U;:  . 

oX-'  '1" 

:-..;  .(.f;.  or:  :r.or6  obCxi  x(..  vilbocv.;  oJ:  bl 


a  to 


.■.::;4, '  lol./'  cl a. 'lox-o  ^  no  .ixoq  ' ..  ox'.x.. 
-I- .  ;.: r'O  w.  .1  .iTX.  1  J  '..ro.'.  ::'..i4.‘  .  ^  .I'.txc. li.Li.:  &.d:.'  :..c..  0,110  ‘.'X' 1.  rp-^v'.i. 

1  .  .'  :■  .  b..'::T:.L  a...r  ;•  :...:r.:;M  .-^rrex;  :  :  xb  :'''W:i'  40  .'nr  irroraoxcoxi  ,  u-xfs 
£:  'OX  to-:r. ,  .  xoj  ;  ;;  i.:j:i  ai:r.oo:.j  t  v/b  '.  b  loi:,  .i:to,c.p i-n  o  toiX 


bo  J 1: 


7 


by  straight  tie-lines  and  that  single  phases  are  represented 
by  solid  portions  within  the  prism  since  they  possess  three 
degrees  of  freedom. 
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3.  The  "Ideal  Solution",  "Fugacity" ,  and  "Activity 

Coefficient"  and  their  Application  to  the  Behavior 
of  "Real"  Systems. 

The  distribution  of  a  component  of  a  system  between 
vapor  and  liquid  is  expressed  by  the  equilibrium  constant  K. 
The  equilibrium  constant  is  defined  as  the  ratio  of  the 
mole  fraction  of  the  component  in  the  vapor  phase  Y,  to 
the  mole  fraction  of  the  component  in  the  liquid  phase  x. 

There  are  several  methods  that  can  be  used  to  evaluate 
K. 

Raoult's  law  and  Dalton’s  law  may  be  combined  to 
calculate  concentration  of  each  component  in  the  vapor  and 
liquid  phase. 

Raoult’s  law  and  Dalton's  law  may  be  expressed  as 
follows : 

p  =  X  .  p*  (Raoult's  law) 
p  =  y  .  P  (Dalton ' s  law) 

where  p  =  partial  pressure  of  the  component  in  the  vapor 
phase  as  defined  by  Dalton  or  Raoul t. 
p*=  vapor  pressure  of  the  component  at  the  temperature 
of  the  system . 

P  =  total  pressure  of  the  system. 
x,y  =  mole  fraction  of  the  component  in  the  liquid  and 
vapor  phases. 
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Combining  the  equations  of  these  two  laws  we  have:  xp*  = 

yP  which  can  be  arranged  to  give:  p*  =  y  =  K  where/ 

P  X 

K  is  the  equilibrium  ratio  at  the  particular  temperature 
and  total  pressure.  Two  factors  that  restrict  the  use  of 
equilibrium  constants  calculated  by  Raoult's  and  Dalton's 
laws  are  immediately  apparent: 

(1)  A  pure  compound  cannot  have  a  vapor  pressure  at 
temperatures  above  its  critical  temperature.  Thus, 
Raoult's  law  is  limited  to  temperatures  less  than  the 
critical  temperature  of  all  of  the  compounds  in  the 
system . 

(2)  Dalton's  partial  pressure  law  assumes  that  each 
component  in  the  vapor  phase  behaves  as  an  ideal  gas. 
For  practical  purposes,  the  ideal  gas  assumption 
places  an  upper  limit  of  pressure  of  50  to  100  PSIA. 

A  better  approximation  of  actual  equilibrium  constants 
is  obtained  by  defining  the  constants  in  terms  of  the 
fugacities  of  the  compounds.  The  fugacity  may  be  looked 
upon  as  a  vapor  pressure  modified  to  represent  correctly 
the  escaping  tendency  of  the  molecules,  from  one  phase  into 
the  other.  By  using  fugacities  to  define  equilibrium 
constants,  the  departure  of  the  vapor  from  the  ideal  gas 
behavior^  and  the  effect  of  the  total  system  pressure  on  the 
vapor  pressure  of  the  components^  are  taken  into  account. 
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In  an  ideal  solution  the  fugacity  of  any  component  in  the 
system  is  equal  to  the  fugacity  of  the  component  in  the 
pure  state  times  its  mole  fraction.  This  is  frequently 
referred  to  as  the  Lewis  and  Randall  rule. 

Thus  if  fjj  and  fy  are  the  fugacities  of  a  component 
in  the  liquid  and  vapor  phases  respectively  and  f^^  and  fy 
in  the  pure  states 


At  equilibrium  conditions  the  fugacities  of  the  components 
are  equal  in  the  liquid  and  vapor  phase,  so  that  (eq.  1) 

— ^  ■■  ~  =  K  (ideal)  .  Ideal  constants  are  limited  by  one 

X  fv 

important  fact: 

At  a  given  temperature  a  compound  can  exist  as  a 
single  phase  at  any  pressure  other  than  the  vapor  pressure, 
so  that  the  fugacity  of  either  the  vapor  or  liquid  (depending 
upon  whether  the  pressure  is  greater  or  less  than  the  vapor 
pressure  of  the  component)  must  be  obtained  by  extrapolation 
of  the  fugacity  pressure  relationship. 

To  obtain  numerical  values  for  fy  and  fj^  to  use  in 
relation  {J)  one  has  to  integrate  graphically  the  relation: 


(eq.  2) 


that  is  integrate  the  specific  volume-pressure  curve  from  the 
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vapor  pressure  p*  to  the  operating  pressure  P  at  the  fixed 
temperature  of  the  mixture. 

As  the  above  mentioned  extrapolation  is  somewhat 
arbitrary  the  farther  the  system  pressure  is  from  the 
vapor  pressure,  the  more  uncertain  the  equilibrium  constant 
becomes.  Usually  the  ideal  equilibrium  constants  are  considered 
undependable  at  pressures  more  than  twice  the  vapor  pressure 
of  the  compound.  For  this  reason,  they  have  limited  use¬ 
fulness  for  calculations  dealing  with  petroleum  reservoir 
conditions . 

In  considering  the  behavior  of  a  "real"  solution,  it 
is  more  convenient  to  use  the  activity,  activity  coefficient 
and  activity  coefficient  ratio  concepts.  Lewis  and  Randall 
defined  the  activity  of  a  component  in  the  vapor  phase  of  a 
system  as  the  ratio  of  the  fugacity  of  the  component  in  the 
vapor  phase  over  the  fugacity  of  the  pure  component  at  the 
temperature  and  total  pressure  of  the  system. and, 
also,  the  activity  of  the  same  component  in  the  liquid  phase, 
as  the  ratio  of  the  fugacity  of  the  component  in  the  liquid 
phase  over  the  fugacity  of  the  pure  component  as  a  liquid  at 
the  system  temperature  and  the  corresponding  vapor  pressure 
of  the  component  in  the  system. 

In  a  mathematical  form  these  two  expressions  are: 


fv 

av  = 


fL 


f-p 


(eq.3) and  a^ 


f* 


(eq.4  ) 
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Dividing  equation  (  3  )  by  the  mole  fraction  in  the 

vapor  phase,  we  obtain  the  activity  coefficient  of  the 

component  in  the  vapor  phase  of  the  system:  'iv  =  —  ■  (eq.  5) 

Y 

Dividing  equation  (  4  )  by  the  mole  fraction  in  the  liquid 
phase,  etc.,  we  obtain  the  activity  coefficient  in  the 
liquid  phase  L  =  (S'!-  6) 

\J  K  / 

The  ratio  of  Q  ^  =|J(eq.7)  is  the  activity  coefficient 

)( ^ 

ratio.  Substituting  values  and 

(aq.  5)  into  (eqf7)  we  obtain: 


•V  


or 


=  K 


xf  * 


X 


f* 


since 


yfp  -  -P 

_  _  f 

at  equilibrium  we  have  f^  =  f^. Apparently,  represents  a 
correction  factor  for  the  deviation  from  the  real  behavior, 
when  the  fugacity  method  is  applied  in  the  determination  of  K. 
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C.  EXPERIMENTAL  TECHNIQUES 

The  lack  of  data  in  phase  Equilibrium  is  primarily 
due  to  the  difficulty  in  the  construction  and  operation 
of  equipment  which  is  sufficiently  adequate  to  yield 
quantitative  phase  data.  Although  equations  are  available 
to  predict  the  composition  of  phases  existing  at  equilibrium, 
they  are  limited  in  scope  by  assumptions  employed  in  their 
derivation . 

Experimental  data  are  also  required  to  determine 
constants  in  the  equations  for  a  particular  system.  Therefore, 
although  an  extensive  experimental  program  for  a  given 
system  is  time-consuming  it  is  the  best  means  now  available 
for  the  determination  of  phase  equilibrium. 

Although  the  aims  of  contemporary  research  in  the 
field  of  phase  behaviour  are  similar,  the  types  of  equipment 
in  use  differ  considerably.  The  ultimate,  purpose  of 
each  however,  is  to  facilitate  the  control  and  measurements 
of  the  primary  phase  variables  such  as,  pressure,  temperature 
and  phase  composition. 

(5) 

Scheeline  and  Gilliland  clearly  state  the  basis  for 
two  common  methods  and  point  out  the  features  of  each. 

i.  THE  CONSTANT  VOLUME  TECHNIQUE 

In  this  method  the  mixture  of  hydro  carbon  is  introduced 
into  a  high-pressure  bomb,  which  can  be  placed  bodily  in  a 
thermostatically  controlled  bath  of  any  desired  temperature. 

The  apparatus  is  equiped  with  a  stirring  or  shaking  device. 
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vapor  and  liquid  sampling  lines^  and  temperature  and  presure 
indicators.  Agitation  is  continued  until  equilibrium 
between  the  phases  is  reached,  at  which  time  samples  of 
liquid  and  vapor  are  with  drawn  and  analyzed.  The 
method  suffers  from  a  number  of  drawbacks  which  tend  to 
complicate  the  apparatus.  It  is  essential  that  two  phases 
be  present  at  equilibrium,  therefore  the  correct  amount 
of  each  component  must  be  added  to  the  bomb  at  the  beginning 
of  the  run  and  some  level  indicating  device  is  necessary. 

As  critical  conditions  are  approached  and  vapor  and  liquid 
properties  become  similar,  it  is  difficult  to  control 
and  measure  the  phase  interface.  In  order  to  prevent 
pressure  reduction  and  corresponding  flashing  of  liquid 
when  samples  are  with^drawn,  it  is  necessary  either  to 
isolate  the  two  phases  in  separate  compartments  o€  the 
bomb  or  to  add  an  inert  substance  (usually  mercury)  as  the 
sample  is  taken,  which  is  equal  in  volume  to  the  hydro  carbon 
removed.  In  order  to  ensure  the  equilibrium,  agitation  must 
be  prolonged  and  thorough,  but  proviaions  must  be  made 
for  preventing  liquid  droplets  from  entering  the  vapor 
sampling  line. 

2.  THE  VARIABLE  VOLUME  TECHNIQUE 

The  second  method  consists  of  obtaining  the  dew  point 
and  bubble  point  curves  of  a  number  of  prepared  samples  of 
known  composition.  This  can  be  done  either  by  observing  a 
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break  in  plot  of  temperature  vs  pressure  for  a  given 
sample  or  by  enclosing  the  sample  in  a  thick-walled  quartz 
tube  in  a  thermostat  and  observing  the  pressure  at  which  the 
first  trace  of  vapor  or  liquid  appears.  By  cross-plotting 
such  dew  and  bubble  point  data,  points  are  found  at  which 
saturated  liquid  and  saturated  vapor  of  different  compositions 
exist  at  the  same  temperature  and  pressure.  Such  points  in¬ 
dicate  phase  equilibrium.  One  advantage  of  the  method  is 
that  if  conveniently  gives  valuable  data  on  specific  volumes, 
however  a  number  of  serious  drawbacks  limit  its  general 
usefulness . 

The  material  used  must  be  carefully  purified  and  free 
in  particular  from  traces  of  fixed  gases,  for  in  the  critical 
region  the  saturation  pressure  is  quite  sensitive  to  small 
amounts  of  air  etc.  Furthermore  an  extremely  large  amount  of 
experimental  work  must  be  done  in  order  to  define  completely 
and  accurately  the  phase  equilibrium  over  all  values  of 
concentration.  The  most  serious  limitation,  of  the  method, 
however,  lies  in  the  fact  that  it  can  be  used  only  for  systems 
of  two  components.  In  more  complex  systems  it  follows 
from  the  phase  rule  that  composition  is  not  a  unique  function 
of  P  and  T  and  that  therefore  dew  and  bubble  point  measurements 
alone  cannot  define  the  composition  of  two  phases  in  equilibrium. 
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3.  The  Experimental  Equipment  and  Technique 

The  apparatus  employed  was  the  one  constructed  by 
J.  A.  Bailey  as  it  had  been  modified  later  by  A.  P.  Lorenzo 
for  his  work.  It  is  a  modification  of  the  apparatus  Professor 
D.  L.  Katz  of  the  University  of  Michigan  used  for  the  phase 
study  of  the  carbon  dioxide-methane  system  (6) . 

It  is  a  variable  volume  apparatus  with  the  principal  unit 
being  the  visual  type  equilibrium  cell  adapted  from  a  Jerguson 
liquid  level  gauge. 

A  similar  cell  was  employed  for  the  measurement  of  the 
desired  amounts  of  the  single  components  added.  Pressure 
changes  were  made  by  the  use  of  a  hand-operated  mercury  pump 
and  electric  coils,  in  the  air  bath  surrounding  the  equilibrium 
cell,  provided  the  necessary  heating  of  the  sample  to  the 
desired  temperature.  An  electronic  relay  in  connection  with 
an  on-off  thermoregulator  held  the  temperature  constantly 
at  the  desired  level. 

Pressures  were  recorded  on  a  Heise  pressure  gauge  and 
temperatures  were  measured  by  the  use  of  copper-constantan 
thermocouples.  The  withdrawl  of  the  gas  phase  after 
equilibrium  had  been  achieved  was  made  at  constant  pressure 
by  mercury  replacement  of  the  gas  leaving  the  cell.  Collection 
of  the  gases  in  a  sample  collector  followed  the  withdrawl 
and  from  there  the  samples  were  taken  to  be  analyzed. 
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After  the  gas  phase  was  entirely  removed  from  the 
equilibrium  cell,  the  pressure  in  the  cell  was  lowered  so  that 
total  evaporation  of  the  remaining  liquid  phase  would  be 
achieved.  The  obtained  vapor  was  then  analyzed  by  employing 
the  same  technique  as  for  the  gas  sample. 

The  Tutweiler  burette  was  used  for  the  determination 
of  hydrogen  sulfide  content,  and  the  Orsat  apparatus  for 
the  methane  content. 

Figure  I  shows  a  schematic  sketch  of  the  apparatus 


used. 
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FIG.  I.  SCHEMATIC  APPARATUS  ASSEMBLY 
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D.  EXPERIMENTAL  RESULTS 

(1)  Phase  equilibrium  data  at  160°F. 

(a)  Methane-hydrogen  sulfide  system 

A  check  on  the  data  reported  in  the  literature  for 
the  methane-hydrogen  sulfide  system  at  1330  PSIA  and  160°F. 
was  made  and  good  agreement  was  achieved  as  it  is  shown 
below: 

P  =  1330  PSIA  T  =  160°F. 

Mole  fraction  Data  of  Reamer 

of  CH^:  Sage  and  Lacey  This  work 

y  .270  .267 

X  .085  .082 

(b)  The  ternary  system 

The  phase  equilibrium  at  160®F.  has  been  described  by 
the  experimental  determination  of  the  dew  point  and  bubble 
point  loci  at  pressures  of  1000  and  1600  PSIA  by  Lorenzo 
and  at  1330  PSIA  in  this  work.  The  experimental  data  at 
1330  PSIA  are  reported  in  Table  I  in  the  appendix.  The 
isobaric-isothermal  phase  diagram  in  Figure  II  illustrates 
the  equilibrium  at  the  pressure  of  1330  PSIA  which  lies 
between  the  critical  pressure  of  the  carbon  dioxide-hydrogen 
sulfide  binary  system  (1220  PSIA)  and  of  the  me thane -hydro gen 
sulfide  binary  system  (1660  PSIA) . 
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Figure  III  illustrates  the  phase  equilibrium  of  the 
system  at  1300  PSIA  and  160°F.  as  it  was  predicted  by 
Lorenzo  using  his  data  at  1000  PSIA  and  1600  PSIA  as  well 
as  published  data  on  the  binary  mixtures. 

(2)  Equilibrium  ratios 

Experimental  equilibrium  ratios  for  each  of  the  three 
components  were  calculated  at  regular  values  of  the 
parameter  C  up  to  0.6  where  the  two  phase  loop  eventually 
terminates.  Table  II  in  the  appendix  represents  the 
experimental  equilibrium  ratios. 

A  substantial  difference  between  the  experimental  K 
values  for  the  three  components  at  the  pressure  of  1330 
PSIA  and  the  K  values,  picked  from  the  pressure  vs  K  plots 
in  Lorenzo *s  work  is  noted.  Theoretical  K  values  were 
also  calculated  employing  equation  (2)  for  each  of  the 
three  components. 

For  the  determination  of  the  vapor  pressure  of  methane 
and  carbon  dioxide,  both  being  over  their  critical  temperature 
at  160®F.,  an  extrapolation  of  their  vapor  pressure  vs 
temperature  plots  had  to  be  made.  The  reaulting  vapor  pressures 
were:  methane:  9300  PSIA  and  carbon  dioxide:  2350  PSIA.  The 
hydrogen  sulfide  being  below  its  critical  temperature  at  160°F., 
exhibits  a  vapor  pressure  of  779  PSIA. 
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FIG.  3  PREDICTED  MOLAR  COMPOSITION  DIAGRAM 
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Pressure-volume  data  were  available  for  the  pure 
components  methane  and  carbon  dioxide  in  the  literature 
(15)  and  (4  )  respectively.  Values  of  theoretical  K  for 
hydrogen  sulfide  at  160°F.  were  obtained  from  generalized 
fugacity  charts  (19)  rather  than  by  extrapolation. 
Theoretical  K  values  at  regular  C  values  calculated  as  above 
are  lasted  in  Table  III  in  the  appendix. 

(3)  Activity  coefficients  ratios 

The  calculated  theoretical  K  values  were  used  to 
determine  the  activity  coefficient  ratios  of  the  components 
in  the  system  at  160®F .  and  pressure  of  1330  PSIA  which  is 
indicative  of  the  deviation  of  the  liquid  and  vapor  phases 
from  the  ideal  behavior.  Table  III  in  the  appendix  shows 
the  values  of  the  calculated  activity  coefficient  ratios  of 
the  components  at  regular  C  values. 

Values  of  the  activity  coefficient  ratios  of  less  than 
unity  indicate  that  the  component  actually  exists  in  higher 
concentration  in  the  liquid  phase  and  lower  concentration  in 
the  vapor  phase  than  the  ideal  coefficient  ratios  predict. 
Activity  coefficient  ratios  greater  than  unity  indicate 
concentration  of  the  component  in  vapor  phase  greater  than 
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E.  CONCLUSIONS 

In  an  effort  to  discuss  the  accuracy  of  the  obtained 
results,  one  should  mention  several  factors  that  played 
a  significant  role  in  their  achievement. 

1.  Temperature  control  at  the  temperature  of  160°F.  was 
a  difficult  matter  because  of  the  nature  of  the  bath 
and  the  difference  intemperature  with  the  surroundings. 

A  temperature  gradient  of  1°F.  was  always  apparent 
between  the  top  and  bottom  parts  of  the  cell. 

2.  Pressure  control  was  rather  simple  with  an  accuracy  of 
±2  PSIA  for  the  high  pressure  range. 

3.  Volume  measurements  should  be  considered  of  an  accuracy 
of  better  than  ±0.5%.  However,  an  uncertainty  in  the 
overall  system  composition  remains  because  of  several 
traps  of  gases  in  the  way  from  the  measuring  cell  to  the 
equilibrium  cell. 

4.  A  serious  cause  of  mistakes  is  due  to  the  difficulty  in 
separating  the  phases  at  equilibrium.  Condensation 

due  to  throttling  during  the  withdrawal  of  the  gas  phase 
and  sudden  pressure  reduction  destroying  the  achieved 
equilibrium  necessitated  repeats  for  several  runs. 
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5.  The  analytical  method  could  reproduce  methane 

concentrations  (Or sat  apparatus)  of  better  than  ±1.0%. 
The  hydrogen  sulfide  titration,  using  the  Tutweiler 
burette,  gave  an  overall  ±1.0%  error,  achieved  only  by 
duplicating  samples. 

Although  an  overall  error  of  ±1%  in  the  phase 
composition  could  cause  a  substantial  uncertainty  in  K 
values,  and  particularly  for  components  in  very  small 
concentrations  in  the  phases,  the  obtained  results  give  a 
very  good  outline  of  the  phase  composition  for  the  ternary 
mixture. 
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PART  II.  VOLUMETRIC  BEHAVIOR  OF  THE  METHANE- 
CARBON  DIOXIDE -HYDROGEN  SULFIDE 
SYSTEM  IN  THE  VAPOR  PHASE  REGION 
AT  TEMPERATURES  OF  160°,  100°  and 


40°F.  AND  PRESSURES  UP  TO  3000  PSIA 
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F.  INTRODUCTION  -  THEORY 

The  primary  purpose  of  the  experimental  program 
carried  out  in  this  main  part  of  the  thesis  was  an  in¬ 
vestigation  into  the  volumetric  behavior  of  the  methane- 
carbon  dioxide-hydrogen  sulfide  system  in  the  single  (vapor) 
phase  region.  Thus,  the  compressibility  factor  of  a  number 
of  mixtures  representative  of  the  single  phase  (vapor)  region 
was  measured  at  temperatures  of  40,  100  and  160®F.  and 
at  pressures  up  to  3000  PSIA. 

The  compressibility  factors  of  several  mixtures 
of  varied  composition  of  the  carbon  dioxide-hydrogen  sulfide 
binary  system  at  160°F.  were  also  investigated. 

Finally,  an  attempt  was  made  to  compare  the  obtained 
compressibility  factors  with  the  ones  calculated  by  using 
various  correlation  methods  for  mixtures  such  as  the 
Amagat's,  Gilliland's  and  Kay's  methods. 

The  basic  equation  expressing  the  PVT  behavior  of 
gases  is  the  ideal  law: 

PV  =  — ^  (eq.a 

Thus  as  expressed  by  equation  (8  ) ,  the  volume  of  a  unit 
weight  of  an  ideal  gas  at  any  equilibrium  temperature  and 
pressure  can  be  evaluated  solely  from  a  knowledge  of  the 
molecular  weight  of  the  gas. 
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In  practice,  it  has  been  found  that  no  gas  obeys  this 
simple  law  over  all  pressures  and  temperatures. 

The  molecules  of  real  (as  distinct  from  ideal)  gases 
occupy  a  certain  definite  volume  that  may  or  may  not  be 
negligible  in  comparison  with  the  total  volume  occupied  by 
the  gas;  furthermore,  they  possess  fields  of  force  resulting 
in  attractive  forces  between  them.  Both  these  conditions 
lead  to  the  fact  that  as  the  molecules  of  a  gas  approach  one 
another  (or  as  the  gas  becomes  denser)  the  gas  departs  more 
and  more  from  the  ideal  state. 

In  a  few  cases,  unusually  large  deviations  from  the 
ideal-gas  laws  occur  even  at  relatively  low  pressures. 

These  are  thought  to  be  due  to  chemical  changes  resulting 
in  association  or  dissociation  of  molecules.  Hydrocarbons, 
in  particular,  deviate  widely  from  the  ideal  gas  law  at  the 
elevated  pressures  and  temperatures  of  petroleum  reservoirs. 

In  order  to  express  a  more  exact  relationship  between 
the  variables  P,  V  and  T  a  correction  factor  must  be  introduced 
into  equation  (8  )  . 

Either  of  two  methods  is  normally  used:  The  first  method 
introduces  a  multiplying  factor  into  the  right  side  of  the 
equation  (  8  )  and  the  second  method  makes  use  of  an  addition 
factor.  The  multiplying  factor  goes  by  several  names;  three 
of  the  most  common  are  "compressibility"  factor,  "deviation" 
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factor  and  " super -expansibility"  factor.  The  cornpressibility 
factor  and  deviation  factor  are  synonymous,  and  the  super¬ 
expansibility  factor  is  the  reciprocal  of  these  two  factors. 
The  addition  factor  is  termed  the  "residual  volume".  Both 
methods  are  accurate  to  the  same  degree,  but  the  first  is 
easier  to  use  and  therefore  the  most  popular.  The  "residual 
volume"  correlation,  however,  is  of  a  wide  use  in  the 
determination  of  other  thermodynamic  properties  of  mixtures 
employing  PVT  data. 

The  compressibility  factors  of  different  gaseous 
components  vary  over  wide  ranges  of  elevated  pressures. 

During  the  1860 ‘s  many  papers  appeared  in  the  literature 
concerning  PVT  behavior  of  individual  gases.  J.  D.  van 
der  Waals  (7)  in  1873  developed  the  theorem  of  corresponding 
states  on  which  the  modern  correlation  of  compressibility 
factors  is  based.  The  theorem  proposes  that  all  gases  will 
exhibit  the  same  behavior  when  viewed  in  terms  of  reduced 
pressure,  reduced  volume,  and  reduced  temperature.  The 
term  "reduced"  means  that  each  variable  is  expressed  as  a 
ratio  of  the  critical  value.  Thus,  all  variables  are 
reduced  to  a  common  datum,  the  value  of  the  critical.  The 
reduced  conditions  of  volume,  temperature  and  pressure  are 
expressed  as: 
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The  theorem  of  corresponding  states  when  applied  to 
gases  having  similar  chemical  structure  (as  the  paraffin 
hydrocarbons) ,  offers  a  correlation  method  with 
satisfactory  accuracy  for  engineering  work.  However,  the 
theorem  of  corresponding  states  is  not  perfect  leading  to 
a  greater  deviation  when  dealing  with  gases  of  different 
molecular  structure. 

The  importance  of  either  2  or  oC  lies  in  the  ability 
to  correlate  relatively  accurately  the  values  for  different 
gases  by  use  of  the  reduced  conditions  and  the  consequent 
elimination  of  use  of  complicated  equations  of  state 
involving  a  multitude  of  constants  or  higher  degree 
equations,  with  virial  coefficients,  such  as  found  in  the 
theoretical  van  der  Waals  and  Beattie-Bridgmann  or  the 
empirical  Maron  and  Turnbull  or  Benedict-Webb-Rubin 
equations . 

Compressibility  of  Mixtures: 

Although  a  relatively  accurate  method  for  predicting 
compressibility  factors  of  the  single  components  is  provided 
through  the  use  of  reduced  conditions  -  compressibility  charts, 
determination  of  the  correct  values  of  compressibility  factors 
for  mixtures  is  somewhat  difficult. 

Two  general  methods  of  dealing  with  gaseous  mixtures 
have  been  proposed.  The  first  assumes  additivity  direct  or 
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modified  of  certain  physical  properties  of  the  mixture  in 
terms  of  those  of  the  components;  the  second  averages  the 
confetants  of  an  equation  of  state  applicable  to  the  pure 
components.  The  familiar  illustrations  of  the  first  method 
are  Dalton's  (11)  and  Amagat^s  ( ll)  laws.  It  has  been 
generally  assumed  that  Amagat's  law  reproduces  the 
experimental  data  better  than  Dalton's. 

The  attempt  to  average  constants  in  an  equation  of 
state  goes  back  at  least  to  van  der  Waals.  The  method  has 
been  most  successfully  developed  by  Beattie  (8) ,  but  is 
complicated  for  engineering  work  and  limited  by  lack  of 
the  necessary  constants  for  most  substances. 

Several  other  methods  have  been  proposed  for  the 
treatment  of  problems  involving  mixtures  but  any  or  all  may 
break  down  in  the  critical  region. 

A  more  extensive  use  of  two  of  them,  namely  Kay's  rule 
and  Gilliland's  method  is  made  for  the  treatment  of  a  variety 
of  mixtures  (9,  10)  . 

Amagat's  method,  and  also  Kay's  and  Gilliland's  method, 
will  be  discussed  in  the  correlation  of  the  experimental 


data. 


:.G  r 


r 


<  ,-..i 


j-ajl:  ■ 


L'l 


:j' 


r 


'■rx.:; 


I 


.C..X  :>  :rU, 


A,' ;  ■ 


J  '  J-X 


■iV.- 


V  G 


G  "J- 


:iGf- 


.‘iUC  G 


:iS. 


Ll- 


...i 


z'-'f..  G-/'/  . 


VGG  XX.  ..  'I  ■  vXi 


f. 


Li  - 


.  £j.  . . 


r.,A' 


G-r-- 


DG".  '.;J " 


^  G 


i.-. 


•i 


j 


31 


O.  EXPERIMENTAL  TECHNIQUE 

The  experimental  procedure  consisted  of  confining 
a  known  mass  of  sample,  of  definite  composition,  over 
mercury  in  the  variable  volume  equilibrium  cell.  The 
sample  was  thermostated,  and  a  known  pressure  was  applied. 

After  equilibrium  was  attained,  the  volume  was  determined 
by  measuring  the  length  of  the  column  of  sample  by  use  of 
the  cathetometer ;  this  length  was  related  to  the  volume  by 
a  prior  calibration  of  the  cell. 

Thus,  each  complete  run  was  separated  into  several 
distinct  portions  which  included  makeup  of  the  desired 
system,  achievement  of  equilibrium  conditions  within  the 
system  and  the  measurement  of  the  variables,  i.e.,  temperature, 
pressure  and  volume: 

1.  Addition  of  Hydrogen  Sulfide: 

With  the  equilibrium  cell  temperature  at  the  desired 
level  the  run  started  with  the  introduction  of  hydrogen 
sulfide  into  it,  after  a  previous  evacuation  of  the  lines 
leading  to  the  equilibrium  cell,  and  consequent  purging  by 
the  hydrogen  sulfide.  After  a  brief  time  interval  to  allow 
the  hydrogen  sulfide  to  reach  thermal  equilibrium,  its  tem¬ 
perature,  pressure  and  volume  were  recorded  and  thus  a 
calculation  of  the  number  of  moles  of  hydrogen  sulfide 
added  was  possible. 
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2 .  Carbon  dioxide  and  methane  had  to  be  determined  individually 
in  the  measuring  cell.  Since  the  number  of  moles  of  carbon 
dioxide  and  methane  was  fixed  by  the  knowledge  of  the  desired 
overall  system  composition  and  the  number  of  moles  of 
hydrogen  sulfide  already  in  the  equilibrium  cell,  the  volumes 

of  each  component  needed  in  the  measuring  cell  could  be 
predetermined  for  the  pressure  at  which  they  were  to  be 
added  and  temperature  of  the  measuring  cell.  The  order  in 
which  methane  and  carbon  dioxide  were  measured  and  introduced 
into  the  equilibrium  cell  was  not  fixed;  the  procedure, 
however,  was  the  same  for  each  component.  After  the  desired 
amount  of  each  of  the  single  components  had  been  added  and 
their  number  of  moles  had  been  calculated,  the  single 
component  was  injected  into  the  equilibrium  cell  by  mercury 
displacement  from  the  measuring  cell. 

3.  The  temperature  in  the  equilibrium  cell  was  measured 

with  3  '  one- junction  copper-constantan  thermocouples  and  a 

potentiometer  reading  to  2.5  microvolts.  A  one- junction 
similar  thermocouple  was  employed  for  the  measuring  cell. 

4.  pressures  were  determined  with  a  0-3000  PSIA  Heise 
gauge  for  the  higher  pressure  and  with  both  the  above 
mentioned  gauge  and  a  0-1000  PSIA  for  the  lower  pressure 
range.  The  gauges  were  calibrated  against  a  dead-weight 
piston  gauge  and  proper  calibration  curves  were  employed. 
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Suitable  corrections  were  made  for  the  various  hydrostatic 
heads  existing  between  the  sample  and  the  gauges. 

5.  A  cathetometer  was  used  to  measure  the  sample  volume 
and  to  determine  the  hydrostatic  heads  which  contributed  to 
the  measured  pressure.  The  reproducibility  of  the  length 
measurements  was  about  ±  .05  mm.,  corresponding  to  a  precision 
in  volume  measurements  of  about  .01%  for  the  largest  volumes 
observed,  and  about  0.1%  for  the  smallest  ones.  However,  the 
encountered  error  should  be  considered  as  being  greater. 

The  measurements  of  temperature,  total  volume  and 
equilibrium  pressure  were  obtained  upon  both  increase  and 
decrease  in  the  total  volume  of  the  system.  The  agreement 
of  the  values  obtained  by  these  two  approaches  to  a  given 
state  indicated  that  the  attainment  of  a^  steady  state'  was 
satisfactory. 

The  study  was  carried  out  in  the  same  apparatus  that 
had  been  used  for  the  phase  equilibrium  measurements  in  the  first 
part  of  this  work  and  no  modifications  of  it  were  attempted. 
However,  a  more  precise  cathetometer  was  used  for  the  measure¬ 
ment  of  volume  and  refrigeration  coils  replaced  the  heating 
coils  for  the  study  at  40®F. 
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H .  EXPERIMENTAL  RESULTS 

1 .  The  Binary  System  at  160°F. 

The  experimental  program  included  the  investigation 

into  the  volumetric  behavior  of  three  binary  mixtures  of 
carbon  dioxide  and  hydrogen  sulfide  at  the  temparature  of 
160®F.  Table  (  IV  )  in  the  appendix  and  Figure  (  4  ) 

represent  the  results  obtained. 

Since  the  maximum  pressure  during  the  investigation 
did  not  exceed  1800  PSIA  the  compressibility  factor 
experimental  line  was  arbitrarily  extended  up  to  3000  PSIA 
(dash  line) . 

The  investigation  of  the  ternary  system  included  the  study 
of  the  volumetric  behavior  of  several  mixtures  at  the  three 
mentioned  temperatures  of  160,  100  and  40°F.  and  pressure 
up  to  3000  PSIA.  An  attempt  was  made  to  keep  the  methane 
to  carbon  dioxide  mole  ratio  at  4:1  and  vary  the  hydrogen 
sulfide  content;  however,  this  rule  has  not  been  followed 
closely  in  the  preparation  of  the  ternary  mixtures  at  the 
temperatuiE  of  40°F.  since  the  single  phase  (vapor)  region  is 
limited  allowing  the  preparation  of  only  two  mixtures  similar 
to  those  prepared  at  the  temperatures  of  100°  and  160°F. 

2 .  The  Ternary  Mixture  at  160°F. 

At  the  temperature  of  160 °F.  the  following  five  mixtures 
were  investigated: 
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Table 


Mixture  No. 

yH2S 

Yc02 

ycH4 

I 

.0392 

.1814 

.7794 

la 

.0426 

.1816 

.7756 

II 

.1252 

.1751 

.6997 

Ila 

.1118 

.1728 

.7154 

III 

.2331 

.1538 

.6131 

V 

.0857 

.3660 

.5483 

VI 

.5256 

.1044 

.3700 

Since  for  iTiixtures  No.  I  and  II  the  total  number 
of  moles  added  to  the  equilibrium  cell  exerted  a  pressure 
of  well  over  700  PSIA  to  start  with,  i.e.,  with  the  sample 
occupying  the  total  volume  of  the  cell,  two  additional 
mixtures  No.  la  and  Ila^  corresponding  respectively  to 
mixtures  No.  I  and  II^  had  to  be  prepared  for  the  range 
of  pressure  between  0  and  700  PSIA. 

To  cover  the  entire  range  of  pressures  in  the  measure¬ 
ments  of  the  compressibility  of  gases,  a  conventional  method 
is  employed  so  that  the  same  mixture  is  kept  in  use  for  both 
high  and  low  pressure  ranges.  It  consists  of  adding  a  relatively 
great  number  of  moles  into  the  equilibrium  cell  and  after 
measurements  have  been  obtained  for  the  high  pressure  range, 
releasing  part  of  the  sample  so  that  the  remaining  gas 
would  permit  measurements  at  lower  pressures. 


:.|fca' 


V,.  .\ 


Ir  - ;.T 


.,1,. 


T 


Cf.:.  cf::.',,/  ^  .::ixvv’  :i-.;xx:''a  xf.xvc  . ^ 

..r.-  ,.c  j:::  .U'.. cx.i.  ..^.w  ^  .  ..aa.’.'’ .[..-xvC'j  .  .Ta  ,  i 

.cd-  Y-^  '  .C/.  i  ■  ■'-  ’  -  a  . ^  -Xu 

3i:..n."^:£  ...xli  d  .  l.'.v -vi  '  --  ■■ 

■„  ^.•ev  c  ■ /a...  X x;'x,  :.o 


xidd-  :-.ox;.-a....£v;Y  .:■.'...  ■-'  3X.tJ'.cL-  .xd'::'  a.'V^ 

:i:x..:.;  .a:  .a  a.su  V.'.' ._.a. 


ucJ  c,..t  a.i7:;xj:ia  ...xxai';  ■  add'  d'a.dJ  ua  "xc  a.-.: 

.Y.--.:  Ye  Y:..  ...lixckx,;..  J  1  exx  xc  I  a/::.  liY-'-- 

x:,  :.':,,:  /X-...  !..jj.Ya  ',.LLajrx  Yr'.  ..tofir  Yo'Xx.  xa,;..: 

^  i.i.1  '::  x/^icu;  xc,  ..■..,;J  c.  iY:.  xod  dfrj,.;  ja.;Yc[<  ixjaci  -'ii  t ,'/....  ..x .XiX 


ar-j  Av' i..  ■:.  ;_.xa  ..;.:,i:'i:;,  xa  •■  -■■Y  ...;  ■  :j '■■'■■  ><i  'oi.^,  a-x:  : 

.  ou/ ■.::./ ad;,  '.a:  . :•  wc,  1  Ya  ad'!:,,  i  x:. X-Jv, .  a  J..k/' xv 


37 


For  the  determination  of  the  number  of  moles  of  the 
remaining  gas,  a  compressibility  factor  taken  from  the 
data  obtained  at  the  high  pressure  region  must  be  employed, 
and  thus,  an  error  introduced  during  the  investigation  at  high 
pressures  would  be  accumulated  in  the  measurements  at  the 
low  pressure  region.  This  reason,  and  the  fact  that  for  a 
length  of  time  an  uncertainty  had  arisen  regarding  the 
validity  of  the  volume  calibration  curves  in  use,  imposed 
the  earlier  mentioned  procedure,  i.e.,  make  up  of  two 
identical  mixtures,  one  for  the  high  pressure  region  and  the 
other  for  the  low  pressures. 

Table  (  v  )  in  the  appendix,  and  Figure  (  5  ) 

represent  the  experimental  compressibility  factors  for  the 
above  mentioned  mixtures.  Mixture  V  for  which  the  methane 
to  carbon  dioxide  mole  ratio  is  other  than  the  stated  4:1, 
consists  an  intermediate  path  in  the  preparation  of  mixture 
la  from  mixture  Ila.  The  maximum  investigation  pressure, 
for  this  mixture,  did  not  exceed  1000  PSIA. 

3 .  The  Ternary  Mixture  at  100°F. 

Establishing  the  new  temperature  of  100®F.  in  the 
equilibrium  cell^  data  were  obtained  for  five  of  the  above 
mentioned  mixtures,  namely,  I,  la,  II,  Ila  and  III. 

Mixture  IV  (yH2S  :  .1989,  yC02  :  .2781;  ycH^  *•  -5230)  was  also 
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investigated,  consisting  a  path  in  the  preparation  of 
Mixture  II  from  mixture  III. 

Table  (VI  )  in  the  appendix  and  figure  (  6  ) 

represent  the  experimental  compressibility  factors  obtained 
at  this  temperature. 

4 .  The  Ternary  Mixture  at  40°F. 

At  this  temperature,  where  the  isobar ic  loops  extend 
considerably  towards  rich  in  carbon  dioxide  and  methane,  and 
poor  in  hydrogen  sulfide  mixtures,  only  two  mixtures,  namely, 
XI  and  XIII  were  studied  having  their  methane  to  carbon 
dioxide  mole  ratio  4:1.  However,  three  other  mixtures  were 
investigated  in  an  attempt  to  cover  the  vapor  phase  region 
considerably. 

The  following  table  indicates  the  composition  of  the 
mixtures  investigated  at  40°F. 

Table 


Mixture  No. 

yH2S 

yco2 

yCH4 

XI 

.2440 

.1575 

.5985 

XII 

.2095 

.2760 

.5145 

XIII 

.1338 

.1762 

.6900 

XV 

.1475 

.2965 

.5560 

XVI 

.1178 

.2365 

.6457 

Mixture  No.^ 

I  -  o 

II  -'a 

III  -XT 

IV 


CHa 


Compressibility  Factor 


Figure  -  6 .  Compressibility 


Pressure  -  pounds  pier  square  inch. 

Factor  Correlation  of  the  Expepmental  Data  for  the  Ternary  Systems  at  100  °F. 
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Mixtures  at  40°F.  were  studied  after  the  investigation 
at  160°F.  and  100°F.  was  entirely  completed,  so  that  new 
mixtures  had  to  be  prepared. 

Table  (  VII  )  in  the  appendix  and  Figure  (  7  ) 

represent  the  experimental  compressibility  factors  obtained 
at  the  temperature  of  40®F. 

Discussion 

A  check  on  the  consistency  of  the  results  obtained 
was  done  by  employing  the  single  components,  methane, 
carbon  dioxide  and  hydrogen  sulfide. 

An  amount  of  a  single  component  was  introduced  into 
either  the  equilibriiam  cell  or  the  measuring  cell  and  after 
the  desired  temperature  was  established,  the  number  of 
moles  of  the  gas  was  calculated  employing  a  compressibility 
factor  taken  from  the  literature.  Then  the  pressure  was 
increased  and  the  new  volume  of  the  gas  in  the  cell  was  recorded. 
The  procedure  was  followed  until  the  maximum  desired  pressure 
was  reached.  Then,  employing  the  as  above  found  number  of 
moles,  the  compressibility  of  the  pure  gas  was  calculated. 

The  obtained  compressibility  factors  were  then  compared  to 
those  reported  in  the  literature. 

For  a  variety  of  temperatures,  in  both  the  equilibrium 
and  the  measuring  cell,  and  a  number  of  approximately  30 
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single  component  runs  and  pressures  up  to  2500  PSIA,  the 
error  encountered  was  less  than  ±  2.5%. 

Another  check  was  made  by  making  up  an  approximately 
50-50  methane  to  carbon  dioxide  binary  system  for  which 
compressibility  data  have  been  published  (13).  ; 

The  agreement  in  the  compressibility  factors  was  better 
than  ±1.0%  and  for  pressures  up  to  3000  PSIA. 
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II.  Correlation  of  the  Experimental  Data 

1.  Amagat's  Method  and  its  Application  to  the  Binary 

and  Ternary  Systems. 

Amagat's  law  leads  to  a  simple  relation: 

Zm  =  ^  +  .  .  .  where  Zm  is 

the  mean  compressibility  factor  of  the  mixture,  Z^,  Zg  .  .  . 

are  compressibility  factors  of  the  pure  components  at  the 
total  pressure  and  temperature  of  the  system  and  Yg  .  .  . 

are  mole  fractions  of  components  A,B  .  .  .  A  disadvantage  of 

the  method  is  that  when  the  conditions  assume  the  liquid  phase 
for  one  or  more  of  the  components  in  the  mixture,  additivity 
of  volumes  cannot  be  assumed  anymore.  Tables  {VIII,  JX,  X)  in 
the  appendix  and  Figures  (8, ’9, 10)  show  a  comparison  of  the 
experimentally  obtained  compressibility  factors  of  several 
ternary  mixtures,  to  those  obtained  by  applying  the 
Amagat  method  of  prediction  to  the  same  mixtures  at  the  same 
conditions.  Compressibility  factors  for  the  pure  components, 
employed  in  these  calculations,  were  taken  from  (13,14)  and  (15,16) 
the  reported  experimental  compressibility  factors  of  the  mixtures , 
in  comparison,  have  been  taken  from  smoothed  experimental 
curves . 

Figure  (  8  )  represents  a  typical  plot  of  percent 

deviation  of  values  of  compressibility  factors  calculated  by 
(dash  lines) 

Amagat's  method Afrom  those  obtained  experimentally  (solid 
line) .  The  four  (dash)  lines  represent  mixtures  No.  I,  II,  III 
and  VI  at  the  temperature  of  160°F. 
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Figure -8,  Experimental  And  Calculated  Compressibility  of  the 
Ternary  Mixtures  atl60°F. 
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A  maximum  discrepancy  from  the  experimental  values  is 
apparent  in  this  figure  for  pressures  in  the  neighborhood  of  1000 
PSIA.  It  is  also  apparent  that  the  greater  the  hydrogen  sulfide 
content  in  the  mixture,  the  greater  the  discrepancy  from  the 
experimental  value  of  the  compressibility  factors. 

Looking  at  Table  (VIII)  we  see  that  the  contribution 
by  the  hydrogen  sulfide  to  the  value  of  the  calculated 
compressibility  factor  by  Amagat ' s  method,  at  the  mentioned 
pressure  region  and  temperature  of  160°F.  and  for  all  of  the 

mixtures  is  negligible  compared  to  that  of  the  carbon  dioxide 

and  methane  content.  This  is  due  to  the  fact  that  the 
conditions  impose.^  the  liquid  phase  for  the  hydrogen  sulfide, 
resulting  thus,  in  non-applicability  of  the  Amagat  method,  as 
it  is  stated  previously.  The  same  behavior  is  seen  in  Figure 
(  9  )  representing  mixtures  I,  II,  III  at  the  temperature 

of  100°F.  Here  the  maximum  discrepancy  occurs  at  a  pressure 
in  the  neighborhood  of  500  PSIA  since  the  equilibrium  pressure 

of  hydrogen  sulfide  is  now  only  394  PSIA. 

Comparing  the  apparent  discrepancy  from  the  experimental 
values  of  each  one  of  mixtures  I,  II,  III  at  the  two  temperatures 
and  pressures  up  to  2000  PSIA  (Figures  8,  9  )  we  see  an  increase 

as  the  temperature  decreases.  For  pressures  of  2000  PSIA  to 
3000  PSIA  the  discrepancy  follows  the  temperature  variation. 
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Figure  -  9.  Experimental  And  Calculated  Compressibility  of  the 
Ternary  Systems  at  100  °F. 
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In  Figure  (  10  ) ,  at  the  temperature  of  40°F.,  the  maximum 

discrepancy  occurs  over  the  same  pressure  region  as  for  the 
mixtures  at  100°F.  However,  for  the  very  low  pressure  region  - 
(pressures  between  0  and  500  PSIA  -  the  discrepancy  is  greater 
at  the  lower  temperature  of  40®F.  The  above  conclusions  can  be 
drawn  by  looking  at  the  discrepancy  curves,  in  Figure  (  10  ) , 
of  mixtures  XI  and  XIII,  which  are  similar  to  mixtures  III  and  II 
at  the  temperature  of  100® F.  and  160® F. 

It  is  worth  noticing  that  mixture  No.  VI  of  the 
highest  hydrogen  sulfide  concentration  at  the  temperature  of 
160®F.  and  for  pressures  between  2000  and  2600  PSIA  exhibits 
a  discrepancy  much  lower  than  any  of  the  poorer  in  hydrogen 
sulfide  mixtures  at  the  same  temperature.  For  pressures 
between  2600  and  3000  PSIA  and  also  1400  to  2000  PSIA  its 
discrepancy  is  lower  than  that  of  mixture  No.  Ill,  reversing 
sign  at  the  pressure  of  2300  PSIA. 

A  similar  behavior  is  noticeable  for  mixture  No.  Ill 
at  the  temperature  of  100® F. 

At  the  temperature  of  40®F.  where  for  both  hydrogen  sulfide 
and  carbon  dioxide  and  for  a  wide  pressure  range,  conditions 
impose  the  liquid  phase,  the  discrepancy  is  much  greater 
following  an  order  of  methane  concentration.  The  higher  the 
methane  concentrations,  the  lower  the  discrepancy  from  the 
experimental  values.  Here  also,  and  for  pressures  over  2000 
PSIA,  the  discrepancy  diminishes  substantially  with  the  deviation 
curves  on  either  side  of  the  e^erimental  datum. 
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Figure  -  10.  Experimental  And  Calculated  Compressibility  of  the 
Ternary  Systems  at  40  °F. 
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2.  Gilliland's  Method  and  its  Application 
to  the  Binary  and  Ternary  Systems. 

One  of  the  classical  approaches  to  the  study  of  the 
behavior  of  gases  is  the  assumption  that  the  total  volume 
occupied  by  the  gas  may  be  considered  as  made  up  of  the 
effective  volume  of  the  molecules  themselves  and  of  the  free 
volume  available  for  molecular  motion,  and  similarly  that  the 
total  pressure  is  the  sum  of  the  actual  pressure  externally 
exerted  upon  the  gas  and  the  so-called  internal  pressure, 
equivalent  to  the  mutual  attraction  of  the  molecules. 

Also,  one  of  the  most  useful  correlations  of  the  PVT 
properties  of  pure  gases  has  been  the  graphical  representation 
of  the  isometrics,  i.e.,  a  graph  of  pressure  vs  temperature 
at  constant  volume.  In  the  range  above  the  critical 
temperature  the  isometrics  are  essentially  straight  lines  for 
a  large  number  of  substances.  Lewis  and  Luke  (  17  )  have 

expressed  the  equation  of  the  isometric  as  P  =  02.  T  -  02 
with  02^  f  02  being  only  functions  of  volume  if  the  isometrics 
are  straight  lines.  The  reciprocal  of  the  slope  of  the 
isometric  (  l/0i  )  is  a  measure  of  the  free  volume  and  the 
negative  intercept  of  the  tangent  to  the  isometric  (  02  ) 

is  generally  defined  as  the  internal  pressure. 

For  mixtures,  Gilliland  presumes  that  the  above 
equation  can  be  written  as: 
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^rriix  -  (^l)mix  “  (^2^ mix  where 

(^l)m  =  Ya  (^i)a  +  Yb(^^i)b  +  •  • 


mix 


where 


(^2 )  m  -  (Ya  J  (^2^  A  +  Yb  J  (02^ B  ^  ^  ^ 


For  binary  systems  the  equation  could  be  written 


as ; 


(eq.9)  Pmix 


The  values  of  P^,  Pg  are  pressures  exerted  by  the  pure  compo¬ 
nents  alone  under  the  existing  temperature  and  molal  volume. 
Values  of  (02)  A'  (02)  b  calculated  from  a  plot  of 

vs  (02) r  for  each  component.  Equation  (  9  )  becomes 

more  complicated  when  more  than  two  components  are  present. 

The  method  is  a  very  tedious  one  and  hence  of  little 
use  at  present.  An  attempt  was  made  to  compare  the 
experimentally  obtained  data  for  the  hydrogen  sulfide-carbon 
dioxide  system  to  values  calculated  by  this  method.  The 
agreement  between  the  calculated  compressibility  for  the 
mixture  and  the  experimental  was  comparable  to  that  obtained 
by  using  Kay's  method  for  the  ternary  system  as  it  will  be 
shown  later.  An  outline  of  the  method  employed  is  as  follows 

From  pressure  vs  compressibility  factor  tables  for  the 
pure  components  carbon  dioxide  (4  )  and  hydrogen  sulfide  (14) 
the  isobaric,  i.e.,  temperature  vs  molal  volumes  at  constant 
pressures,  lines  are  plotted.  Pairs  of  values  of  temperature 
and  pressure  for  equal  volumes  are  picked  on  the  isobarics 
which  are  then  plotted  on  a  temperature  vs  pressure  graph 
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forming  the  so-called  isometrics.  The  isometrics  for  carbon 
dioxide  are  very  well  straight  lines,  the  isometrics  for 
hydrogen  sulfide  however  are  lines  with  a  slight  curvature. 
Extrapolation  of  the  carbon  dioxide  straight  isometrics 
until  they  cut  the  pressure  ordinate  then  follows.  Extrapolating 
the  straight  part  of  the  hydrogen  sulfide  isometrics,  we  are 
also  able  to  obtain  the  intercepts  with  the  pressure  ordinate. 
Since  the  negative  values  of  this  intercept  represent 
value  ^2  the  isometrics  we  then  tabulate  V  (isometrics) 
vs  02  carbon  dioxide  and  for  hydrogen  sulfide.  Dividing 

V  by  the  critical  volume  of  each  of  the  components  and  02 
by  the  critical  pressure  of  each  of  the  components  we  finally 
evaluate  V^.  and  {02^^  r  each  of  the  components.  The  vs 

(02)2:  correlation  is  shown  for  hydrogen  sulfide  and  carbon 
dioxide  in  Figure  (11)  in  the  appendix.  The  method  then 
follows  a  trial  and  error  procedure  which  is  thoroughly  developed 
in  the  appendix  (page  84)  and  for  the  compressibility  factor 
of  carbon  dioxide-hydrogen  sulfide  mixture  No.  2  at  1000  PSIA 
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3.  Kay's  Method  and  its  Application  to  the  Binary 

and  Ternary  Systems. 

As  it  was  stated  previously,  the  modern  form  of  the 
theorem  of  corresponding  states  covers  nearly  all  pure  com¬ 
pounds  so  that  the  calculation  of  their  compressibility  at 
any  condition  is  not  a  difficult  matter  anymore.  If  the  data 
on  complex  mixtures  are  plotted  in  a  similar  manner  by 
employing  the  actual  critical  pressure  and  temperature  of 
the  mixture  for  the  calculation  of  the  reduced  properties, 
the  resulting  curves  deviate  considerably  from  those  of 
pure  compounds,  particularly  in  the  critical  region. 

W.  B.  Kay  stated  that  this  is  due  to  the  improper  selection 
of  the  critical  point. 

Obviously,  if  the  data  in  the  superheated  region  are  to 
be  compared  with  those  for  a  pure  compound  by  the  theorem  of 
corresponding  states,  the  logical  point  for  calculating  the 
reduced  quantities  is  not  the  actual  critical  point  but  a 
point  which  is  the  critical  point  of  the  pure  compound  whose 
PVT  relations  in  the  superheated  region  are  identical  with 
those  for  the  mixture.  If,  then,  the  pressure  and  temperature 
of  this  point  are  known,  the  PVT  relations  in  the  superheated 
region  can  be  calculated  exactly  as  they  would  be  for  a  pure 
compound,  regardless  of  the  boiling  range  of  the  mixture.  This 
point  may  be  called  the  "pseudocritical  point"  of  the  mixture 
since  it  is  analogous  to  the  critical  point  of  a  pure  compound. 
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The  difference  between  this  point  and  the  actual  critical 
point  is  dependent  on  the  width  of  the  boiling  range  of  the 
mixture.  As  the  boiling  range  is  diminished,  the  critical 
point  approaches  the  pseudocritical  and  the  two  points 
coincide  and  are  identical  in  the  case  of  a  pure  compound. 

Since  the  pseudocritical  point  of  a  mixture  has  no 
real  existence,  it  cannot  be  measured.  However,  it  may  be 
estimated  graphically  from  the  PVT  data  on  the  mixture. 

A  thorough  investigation  on  hydrocarbon  mixtures  made 
by  W.  B.  Kay  ( 9  )  yields  the  following  simple  relation  for 
the  pseudocritical  pressure  and  temperature  of  the  mixture: 

(Tc)  mix  =  TcaYA  +  TcbY„  +  (  eq.lO  ) 

a 

(Pc)  mix  =  Pca^A  +  PCgYg  +  (  eq.ll) 

where  Tc^,  TCg  .  .  .  are  critical  temperatures  of  pure  A,  B  . 

.  .  and  Pc^/  Pcg  .  .  .  are  critical  pressures  of  A,  B  .  .  . 

This  method  has  the  advantage  of  using  the  same  methods 
of  solution  of  problems  involving  pure  gases  and  vapors. 

Another  inherent  advantage  of  Kay's  method  is  the  ability  to 
handle  mixtures  where  the  Amagat  rule  would  find  that  the 
conditions  imposed  the  liquid  state  on  one  or  more 
components. 

In  the  original  communication  of  Kay  (9  ) ,  it  is 

remarked  that  for  mixtures  whose  constituents  differ  greatly  in 
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molecular  weight  as  well  as  chemical  nature,  values  of  the 
pseudocritical  pressure  and  temperature,  calculated  as  above 
(eq.lDil)  /  are  likely  to  be  in  error  by  an  amount  greater  than 
that  allowable  in  most  engineering  calculations. 

Figures  (8 , 9  ,  ID)  and  tables  (XT  ,X[E,XIl^in  the  appendix  show 
the  calculated  compressibility  factors  for  most  of  the  studied 
ternary  mixtures.  In  broken  lines  and  in  Figure  (  8  )  the 
discrepancy  between  calculated  and  experimental  values  for 
mixtures  No.  I,  II,  III,  VI  is  shown  at  the  temperature  of 
160°F. 

In  general,  we  see  a  better  approximation  achieved  by 
using  Kay’s  method  over  the  entire  range  of  pressure  than 
Amagat's  at  this  temperature.  However,  although  the  error  is 
much  more  consistent  and  much  lower,  it  exceeds  the  allowable 
for  many  engineering  calculations  proving  the  time-consuming 
experimental  investigations  absolutely  necessary. 

On  Figure  (  9  ) ,  representing  mixtures  No.  I,  II  and 

III  at  the  temperature  of  100°F.,  we  observe  nearly  the  same 
behavior  of  Kay's  method.  However,  an  increasing  discrepancy 
for  mixture  No.  Ill  (highest  hydrogen  sulfide  concentration) 
is  apparent  at  this  temperature  and  for  the  range  of  higher 
pressures,  i.e.,  over  1500  PSIA. 

This  phenomenon  is  more  pronounced  in  Figure  (  10 ) 
representing  ternary  mixtures  at  the  temperature  of  40°F.  At 
this  temperature  and  pressure  over  1000  PSIA  the  discrepancy 
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from  the  experimental  values  of  the  compressibility 
factors  becomes  much  greater  than  that  obtained  using  the 
Amagat's  correlation  method.  With  only  one  exception  (mixture 
No.  1)  the  discrepancy  increases  as  the  methane  content 
decreases  or  as  the  sum  of  hydrogen  sulfide  plus  carbon 
dioxide  increases. 

At  this  temperature,  and  pressures  over  1000  PSIA,  since 
both  hydrogen  sulfide  and  carbon  dioxide  are  in  the  liquid 
phase,  the  systems  depart  greatly  from  ideality. 

This  reason,  and  the  fact  that  both  acidic  constituents 
of  the  ternary  mixtures,  i.e.,  carbon  dioxide  and  hydrogen 
sulfide  differ  considerably  in  molecular  weight  and  chemical 
structure  from  the  main  hydrocarbon  component  of  the  mixture, 
i.e.,  methane,  may  well  explain  the  greater  discrepancy 


observed. 
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4.  Residual  Volume  Correlation 

Asv>as  stated  in  a  previous  chapter,  the  compressibility 
factor,  when  plotted  against  reduced  conditions,  provides  a 
very  satisfactory  correlation  of  PVT  data.  Vapor  pressure 
volumes  of  pure  compounds  estimated  from  this  type  of  plot 
are  of  sufficient  accuracy  for  most  engineering  design  cal¬ 
culations  . 

Kay's  method,  also,  of  the  "pseudoreduced"  conditions, 
provides  a  satisfactory  means  of  correlating  compressibility 
factors  of  mixtures  in  a  similar  manner. 

For  thermodynamic  calculations,  however,  these  correlations 
are  not  accurate  enough  and,  in  addition,  are  not  well 
adapted  to  graphical  methods.  Neither  is  any  of  the  usual 
equations  of  state  sufficiently  accurate  over  wide  ranges  to 
permit  accurate  determination  of  first  or  higher  differential 
coefficients . 

In  the  calculation  of  S,  H,  Cp,  |j,,  etc.,  when  p  and  T 

,  ^V 

are  the  independent  variables,  the  coefficients  (  )p 

and  (  must  be  evaluated.  Deming  and  Shupe  (  18  )  used 

among  other  residual  quantities  the  "residual"  volume  o(  = 

RT  -  V  from  which  several  equations  can  be  obtained  by 
P 

differentiation,  yielding  values  for  the  above  differential 


coefficients . 
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An  inherent  feature  of  the  residual  volume  correlation 

is  its  ability  in  describing  very  accurately  the  behavior 

of  a  gas  at  relatively  low  pressures. 

12,  13  (XV, XVI 

Figures  (  &  14  )  and  Tables  &  XVII)'  in  the  appendix 
represent  such  a  correlation  for  several  of  the  experimentally 
studied  ternary  mixtures.  The  experimental  compressibility 
factors  employed  in  the  calculation  of  the  residual  volume  , 
were  obtained  from  smoothed  experimental  curves. 
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Residual  Volume,  or  cubic  feet  per  pound  mole. 
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Pressure  -  pounds  per  square  Inch. 

Figure- 12.  Residual  Volume  Correlation  of  the  Experimental-  Data 
of  the  Ternary  Systems  at  160  ®F. 


Residual  Volume,  or  —  cubic  feet  per  pound  mole. 


60 


250  500  750  1000  1250  1500  1750  2000  2250  2500  2750  3000 


Pressure  -  pounds  per  square  inch. 

Figure  -  13.  Residual  Volume  Correlatiom  o*f  the  Experimental  Data 
of  the  Ternary  Systems  at  100  °F. 
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Pressure  -  pounds  per  square  inch 

Figure -14.  Residual  Volume  Correlation  of  the  Experimental  Data 
cjf  the  Ternary  Systems  at  40  °F. 
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I.  CONCLUSIONS 

Following  the  discussion  of  the  results  and  the 
comparison  of  the  experimental  compressibility  factors  to 
those  calculated  by  the  mentioned  methods  a  summary  of  arising 
conclusions  would  have  been  as  follows: 

1.  The  temperature  during  the  investigation  plays  an 
important  role  in  the  question  whether  the  compressibility 
of  the  mixtures  can  be  described  by  means  of  methods  of 
prediction.  Kay's  method  could  be  applicable  with  rather 
satisfactory  results  for  some  engineering  problems  for 
the  temperatures  of  160°  and  100°F .  over  the  entire  range 
of  pressure  studied.  Methods  based  on  additivity  concepts 
are  not  valid  even  for  the  lower  range  of  pressure.  At  the 
temperature  of  40°F.  there  is  no  adequate  method  of  prediction. 

2.  The  pressure,  also,  has  an  important  influence  on  the 
discrepancy  between  the  calculated  compressibilities  and  the 
experimental.  At  higher  pressures  the  discrepancy  is  greater, 
save  the  case  when  Amagat's  method  is  applied,  and  this 
phenomenon  is  more  pronounced  as  the  temperature  decreases. 
Amagat's  method  does  not  follow  the  rule  exhibiting  a  maximum 
discrepancy  over  the  range  of  pressures  for  which  at  the 
temperature  of  the  investigation  one  or  more  of  the  pure 
components  are  assumed  as  being  in  their  two  phase  region 

or  in  the  liquid  phase. 
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3.  The  composition  of  the  mixtures,  also,  is  very 

significant..  Over  the  entire  range  of  pressure 
and  for  the  three  temperatures  of  the  investigation 
it  can  be  said  that  the  discrepancy  increases  as  the 
methane  content  in  the  system  decreases. 

The  limited  number  of  runs,  for  which  an  appreciable 
variation  of  carbon  dioxide  and  hydrogen  sulfide  content 
was  achieved  in  this  work  is  not  sufficient  to  permit 
a  discussion.  A  future  investigation  describing  the 
individual  role  of  hydrogen  sulfide  and  carbon  dioxide 
in  the  compressibility  of  the  ternary  system  would  be 
of  some  significance. 
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APPENDIX 

1.  TABULATED  EXPERIMENTAL  AND  CALCULATED  DATA 

2.  SAMPLE  CALCULATIONS 
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Table  I 

Phase  Composition  Data 
P  =  1330  PSIA  T  =  160°F. 
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Table  II 

Equilibrium  Ratios  at  Regular  Values  of  the  Parameter  C 

P  =  1330  PSIA  T  =  160°F. 
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Table  III 

1.  Experimental  and  Theoretical  Equilibrium  Ratios 

2.  Calculated  Activity  Coefficient  Ratios 
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Table  No.  XV.  Residual  Volume  Correlation  of  the  Ternary  System  at  160°F. 
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Table  No. XVI.  Residual  Volume  Correlation  of  the  Ternary  System  at  100®F. 
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Table  No.  XVII.  Residual  Volume  Correlation  of  the  Ternary  System  at  40°F. 
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Sample  Calculation  -  Part  1 
Run  A  -  8 

Equilibrium  condition  :  T  =  160  Degrees 
Fahrenheit;  P  =  1330  pounds  per  square  inch  absolute 
The  desired  total  composition  in  the  two  phase 
region  was: 

H2S  =  .72  mole  fraction 

CH^  =  .13  mole  fraction 

CO2  -  .15  mole  fraction 


(1)  Addition  of  H2S  to  equilibrium  cell 
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400  psia. 

161  Degrees  Fahrenheit  =  621  Degrees  Rankin 
71.0  c.c. 

0.851 

670  - 

irr  .  °R  .  gm.  mole 

400  X  71.0 

zRT  0.851  X  670  x  621  -0801  gm.  moles. 


(2)  Addition  of  CH4  to  measuring  cell 
(a)  Moles  of  CH4  required: 

n  =  X  .0801  =  .0145  gm.  moles. 

CH4  was  to  be  added  to  the  measuring  cell  at  the  pressure 
of  200  psia  and  at  a  cell  temperature  of  72  degrees  F. 

At  these  conditions  the  compressibility  factor  was  0.975. 
The  volume  required  in  the  measuring  cell  was: 
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nzRT 

P 


0.0145  X  0.975  x  670  x  532 

200 


25.15  cc 


(b)  Actual  quantity  of  CH4  added: 

P  =  200  psia. 

T  =  72  Degrees  F.  =  532  Degrees  R 

V  =  28.2  c.c 


n  = 


PV 


200  X  28.2 


zRT  670  X  0.975  x  532 

(3)  Addition  of  CO2  to  measuring  cell 
(a)  Moles  of  CO2  required: 


=  .0162  gm.  moles 


n 


n  +  .0801  +  X)162 


=  .15 


n  =  .0137  gm.  moles. 

CO2  was  to  be  added  to  the  measuring  cell  at  a 
pressure  of  100  psia.  and  at  a  cell  temperature  of 
72  Degrees  F.  At  these  conditions  the  compressibility 
factor  of  CO2  was  0.958.  The  volume  required  in  the 
measuring  cell,  then,  was: 


V 


nzRT 


.0137  X  0.958  X  670  x  532 

100 


=  46.8  c.c 


(b)  Actual  quantity  of  CO2  added: 

P  =  100  psia. 

T  =  72  Degrees  F.  =  532  Degrees  R 

V  =  60.0  c.c 


n  =  -PV 


100  X  60.0 


zRT 


0.958  X  670  X  532 


=  .0176  gm.  moles 
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System 

Composition : 

Moles 

Mole  Fraction 

H2S  : 

.0801 

.7032 

CH4  : 

.0162 

.1423 

CO2  : 

.0176 

.1545 

(5)  Composition  of  Phases  at  Equilibrium  as  Determined 
by  Analysis: 


H2S 

CH4 

CO2 


Gas  Phase 
.699 
.1395 
.1615 


Liquid  Phase 
.770 
.1075 
.1225 


L>i-^ 


■V.'  I 


C 


CX<.'.  « 


84 


Sample  Calculation  of  Predicting  the  Compressibility  Factor 
by  the  Gilliland's  Method 

Mixture  chosen:  No.  2  binary  of  carbon  dioxide- 
hydrogen  sulfide. 

Temperature  =  160®F.,  Pressure  1000  psia. 

Then,  experimental  compressibility  factor:  . 735 . 

We  assume  a  number  of  compressibility  factors  including 
the  experimentally  found. 

Choose:  z  2_  =  .400 

Z2  =  .600 

Z3  =  .800 

We  then  calculate  the  molal  volume  occupied  by  the  gas 
mixture  of  the  chosen  compressibility  factors  at  the 
above  conditions. 

zt  =  0.4  V-i  =  — ^  ^  Q  ■ .  =  .1662  lit/gmole. 

Z2  =  0.6  V2  =  0^4  ""  -2493  lit/gmole. 

Z3  =  0.8  V3  =  2Vj_  =  .3324  lit/gmole. 

Step  1 

We  then  calculate  the  pressure  exerted  by  each  of  the 
components  of  the  mixture  separately  if  they  were  occupying 
the  above  found  molal  volumes  at  the  temperature  of  160°F. 
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General  equation  to  be  used: 

=  z  — where 
PcV 

R  =  .082  T  =  (160”F.)  or  344°K  and 

=88.9  atm.  Pc''*^2  =  73.0  atm. 

The  procedure  to  obtain  required  a  trial  and 

error  solution  using  the  reduced  conditions  -  compressibility 
charts  (  12  )  for  the  previously  found  values  of  V.  The 

obtained  values  of  P^-  for  each  of  the  components  at  the 
three  values  of  V  are  shown  in  the  following  table: 


Table  a. 


I  CO2 


082  X  344 

73.0  V 


.3864  _z 
V 


-e-  /  T. 


344 

304.1 


1.13 


II  H2S  :  P^ 


.082  X  344  2 

88.9  V 


.3173 


344 

373.4 


.92 


I 


II 


Vi  = 

.1662 

Pri  = 

2 . 32 , 

=  1.35 

V2  = 

.2493 

Pr2  = 

1.55, Z2 

=  1.07 

V3  = 

.3324 

Pr3  ,= 

1.16, Z3 

=  .88 

Vi  = 

.1662 

Pri  = 

1.91, 

=  .70 

V2  = 

.2493 

Pr2  = 

1.27,  Z2 

=  .66 

V3  = 

.3324 

Pr3  = 

•95,  Z3 

=  ,63 

and  thus 


we  obtain  Table  b 
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Table  b. 

V;L  =  .1662 

V2  =  .2493 

V3  =  .3324 

C02 

:  P  =  EV  98.6 

78.1 

64.2 

H2S 

:  P  =  Pr-  Pc  62,2 

58.7 

56.0 

step  2 

From  Figure  (  11  )  ,  representing  Vr  vs(^2We  values 


of 

04 

3-  corresponding  to 

the  ( 

chosen 

values  of 

CO2 

Vrl 

Vl 

.1662 
.0957  ~ 

1.74 

*^2rl 

'  =  2.17 

Vr2 

= 

Vc 

.2493 
.0957  “ 

2.61 

S^2r2 

=  1.05 

Vr3 

V3 

 .3324  

3.47 

i^2r3 

=  .63 

v„ 

.0957 

H2S 

1 - 1 

u 

> 

Vl 

 ol662 

1.70 

S2*2rl 

=  1.280 

Vc 

,0977 

Vr2 

V2 

 .2493 

2.55 

S^2r2 

=  0,300 

Vc 

.0977 

Vr3 

V3 

 .3324 

3.40 

^2x3 

=  0.105 

Vc 

.0977 

So,  we  finally  obtain  values  of  02  each  of  the  components 

for  the  chosen  compressibility  factor S- 


Vl 

=  .1662 

V2  =  .2493 

• 

II 

ro 

> 

CO2  : 

^2 

=  S^2r  Pc 

158.4 

76.7 

46.0 

H2S  : 

02 

=  ^2r  Pc 

114.0 

26.7 

9.4 

3324 


i-og»rithmic.  2  X  2  < >ci.  Ppduced  Volume 
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Reduce  Internal  Pressure  -  (j)2r 

-  (j)^^  Relationship  for  HeS  and  CO2 


Figure  - 11 
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Substituting  the  values  of  P  and  0  2  found 
above  in  the  equation  for  the  binary  mixtures: 

^mix  =  ^002^002  +  Yh2S^H2S  +  yC02^H2S  ^2C02  ~  ) 

we  get:V^  =  .1662  or  =  0-4 

Pmix  =  (.2)x(62.2)  +  (.8)x(98.6)  +  (.2)x(.8)x(  V  (114)  -  'iOSsTI))^ 
=  12.54  +  78.68  +  (.16)X(10.68  -  12.55)^ 

=  12.54  +  78.68  +  ( . 16)x  (-1 .87) ^  =  12.54  +  78.68 
+  0.56  =  91.78  atm. 

V2  =  .2493  or  Z2  =  0.6 

=  (.2)x(58.7)  +  (.8)x(78.1)  +  (.16)x(  \l  26 .7  - 

=  11.74  +  62.48  +  (.16)x(5.16  -  8.80)^ 

=  11.74  +  62.48  +  (.16)x(-3.64)^  =  11.74  +  62.48 
+  0.58  =  74.80  atm. 

V3  =  .3324  or  Z3  =  0.8 

=  (.2)x(56)  +  (.8)x(64.2)  +  (.16)x(  vjl.4  -  1/4O)  ^ 

=  11.2  +  51.36  +  (.16)x(3.06  -  6.78)2 
=  11.2  +  51.36  +  (.16)x(-3.72) ^  =  11.2  +  51.36 
+  .595  =  63.16  atm. 

Plotting  the  assumed  compressibility  factors  vs  the 
as  above  calculated  Pitiix  finally  find  the  compressibility 
factor  at  1000  PSIA  or  68.0  atm.  for  the  system  to  be 
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